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Modification of the classical nucleation theory based on molecular simulation data for 
surface tension, critical nucleus size, and nucleation rate 

Martin Horsch, Jadran Vrabec0 and Hans Hasse 
Institute of Thermodynamics and Thermal Process Engineering, 
Universitdt Stuttgart, Pfaffenwaldring 9, 70569 Stuttgart, Germany 

Nucleation in supersaturated vapor is investigated with two series of molecular dynamics sim- 
ulations in the canonical ensemble. The applied methods are; (a) analysis of critical nuclei at 
moderate supersaturations by simulating equilibria of single droplets with surrounding vapors in 
small systems; (b) simulation of homogeneous nucleation during condensation with large systems 
containing 10^ - 10^ particles for calculating the nucleation rate of vapors at high supersaturations. 
For the Lennard- Jones fluid, truncated and shifted at 2.5 times the size parameter, it is shown that 
the classical nucleation theory underestimates both the nucleation rate and the size of the critical 
nucleus. A surface property corrected modiflcation of this theory is proposed to consistently cover 
data on the surface tension of the curved interface, the critical nucleus size, and the nucleation rate. 

PACS numbers: 64.60.qe, 68.03.Cd, 82.60.Nh 



INTRODUCTION 

Homogeneous nucleation during condensation of super- 
saturated vapors is a well-studied topic; however, it is not 
yet fully understood despite its general importance. The 
most widespread modelingapproach is still the classical 
nucleation theory (CNT) CNT is an accept- 

able approximation for some simple fluids but may yield 
huge deviations compared to experimental data in other 
cases . An important source of error is the assumption 
that the emerging liquid has the same thermodynamic 
properties as the bulk liquid phase @. Condensation 
processes of practical interest, e.g. in atmospheric sci- 
ence, are usually heterogeneous or ion-induced and have 
a more complex mechanism of nucleation Q. However, 
to adequately describe such processes a thorough under- 
standing of the homogeneous case is a prerequisite. 

Experimental methods for studying homogeneous nu- 
cleation face considerable challenges: experimentally, a 
homogeneous system without walls or other irregularities 
can at best be approximated, a difficulty that is absent in 
molecular simulation. Furthermore, the experimentally 
accessible range of the nucleation rate J is limited to com- 
paratively slow processes that are relatively far from the 
spinodal Q. Experimental data on the critical nucleus 
size have only recently become available In molecu- 
lar simulation, homogeneous nucleation can straightfor- 
wardly be studied by a direct approach where a super- 
saturated vapor is observed for some time interval, the 
emerging nuclei are counted, and their size is evaluated 
lol 11 1 . Due to the limitations in computational power, 
accessible system size and time interval are limited. Thus 
the direct approach can currently only be applied to va- 
pors at high supersaturations, where nucleation occurs 
within nanoseconds. 

For systems at lower supersaturation it is necessary 



to follow other, more indirect approaches, e.g. by simu- 
lating other ensembles or related systems instead of nu- 
cleation in the supersaturated phase itself, which occurs 
too slowly. Key quantities determined from such indirect 
simulations are the size of the critical nucleus n* and its 
Gibbs energy of formation AG*, and methods based on 
transition path sampling also permit a study of kinetic 
aspects [12|. Both molecular dynamics (MD) simulation 
with inserted nuclei in noncquilibrium with the surround- 
ing supersaturated phase [13] or based on transitionpath 
sampling [l^ and Monte Carlo (MC) simulation [l^ 
were used for such purposes in the past. 

In the present work, both the direct and an indirect 
simulation approach were applied to validate two ver- 
sions of CNT and to develop a new surface property cor- 
rected (SPC) modification of CNT. A simple model fluid 
was chosen, where the intermolecular interactions are de- 
scribed by the Lennard-Jones (LJ) potential, truncated 
and shifted at an intermolecular distance r = 2.5a [l^ . 
The small cutoff radius leads to relatively fast simula- 
tions and avoids long-range corrections that are hard to 
estimate for inhomogeneous systems [17[ . The truncated 
and shifted LJ potential (LJTS) defines an important and 
well studied model fluid that can be used to describe no- 
ble gases and methane very accurately. A considerable 
amount of thermodynamic data is available for it and 
in particular, the dependence of the surface tension on 
curvature has been quantified [3]. 



CLASSICAL NUCLEATION THEORY 

To describe homogeneous nucleation during condensa- 
tion, a supersaturated vapor in a volume V at the tem- 
perature T and a pressure p which is larger than the 
saturated vapor pressure Ps is considered. The quotient 
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S = p/ps is called the supersaturation of the vapor (with 
respect to pressure). Starting from a homogeneous va- 
por with S > 1, nanoscopic droplets begin to form after 
some induction time as dispersed nuclei of the emerging 
liquid phase. They assume a specific size distribution, 
and the critical nucleus size n* is the number of particles 
n where the Gibbs energy of nucleus formation AGn has 
its maximal value AG* [21 ■ 

The size of the critical nucleus and its energy of for- 
mation were discussed by Gibbs [l| from a theoretical 
standpoint. Given that the number of nuclei with n par- 
ticles is usually determined in CNT by applying a factor 
of exp(— AG,i/fcsT) to the number of monomers, an in- 
ternally consistent approach [l9| leads to the expression 



AG„ 
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Here, C„ is the surface free energy of a nucleus with n 
particles; /ig and /i are the chemical potentials of the sat- 
urated and the supersaturated vapor. In expression ([1]) 
a negative volume contribution competes with a positive 
surface contribution. The difference between the chem- 
ical potentials can be determined from an integral over 
the pressure along the isotherm of the metastable vapor 



P dp 
= P ' 



(2) 



Volmcr and Weber 
by 



approximated the nucleation rate 



J ^Cexp{~AG*/kBT). 
The preexponential coefficient is Q 
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where A* represents the surface area of a critical nucleus, 
A^i the number of monomers constituting the vapor. 
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is the Zel'dovich factor, and m is the mass of a particle. 
Furthermore, Af = b'^ / {b'^ + q^) is the thermal nonac- 
commodation factor which is calculated from 'the energy 
released on addition of a monomer' to a critical nucleus 
'above that needed to maintain the existing temperature' 
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and the kinetic energy variance 

62 = (c, -hfcB/2)fcBT2, 



(6) 
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where Ah^ is the bulk enthalpy of vaporization and is 
the isochoric heat capacity of the vapor [3]. 

CNT is based on the capillarity approximation: the 
density of a nucleus is assumed to be the bulk saturated 
liquid density pi and its surface tension_7„ to be the sur- 
face tension of the planar interface 700 [J, S, B 0] ■ Nuclei 
are assumed to be spherical, thus for one containing n 
particles the surface area is An = [Q-jTin/ piif'^'^ . 

The surface free energy is related to the surface ten- 
sion by 7 = {dC,/ dA)p_T ■ The capillarity approximation 
implies C„ = ^00 An- An analysis of experimental results 
presented by Fenelonov et al. seems to suggest that 
for some fluids, the surface tension of small nuclei de- 
viates from 7oo only by factors between 0.92 and 1.14. 
Their line of argument is based on the so-called 'first 
fundamental nucleation theorem' 21 1 
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according to which the size of the critical nucleus is ob- 
tained from the supersaturation dependence of the nu- 
cleation rate at constant temperature. This value of n* 
is then inserted into the Kelvin equation 
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a corollary of standard CNT, to obtain the surface ten- 
sion of the critical nucleus. However, the Kelvin equation 
does not take any dependence of 7„ on n into account. 
Hence, it is inconsistent to use this equation for quantify- 
ing precisely this size dependence. The nucleation theo- 
rem as given in Eq. ([8]) assumes dp = kBTdlnS , which is 
a bad approximation at high temperatures, in particular 
near the spinodal line. Furthermore, it neglects the de- 
pendence of the preexponential coefficient C from Eqs. 
([3]) and (|4]) on n*, although C is actually proportional 
to both p and the surface area of the critical nucleus (cf. 
Schmelzer [2^ for valid forms of the nucleation theorem). 
From theoretical considerations |23|, |2J| and simula- 



tions @, [1^ it can be inferred that the surface tension for 
interfaces with a high curvature is actually much lower. 
An approximation of the size dependence of the surface 
tension was given by Tolman [23| 



In = 
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where 5 is called the Tolman length and R is the radius of 
the nucleus. Laaksonen, Ford, and Kulmala (LFK) [2^ 
also proposed a size dependent specific surface energy 



Cn/A„ =7oo(l + ain ^'^ + a^n ^Z^), 
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where a-^ and are determined from thermal properties. 
This modification leads to predictions which were found 
to agree better with simulation data of Tanaka et al. 11 1 
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than standard CNT. In the LFK model, the curvature 
effect is covered by the single parameter ai, since AG„ 
does not depend on a^- 



INDIRECT APPROACH: CRITICAL NUCLEI 
FROM MD SIMULATIONS OF EQUILIBRIA 

Phase coexistence methods are an established ap- 
proach for obtaining equilibrium data from molecular 
simulation [26| . In the present indirect simulation ap- 
proach, a single nucleus in equilibrium with a supersatu- 
rated vapor was studied in the canonical ensemble. For 
such simulations it is crucial to choose the relation of the 
number of particles in the nucleus to the total number of 
particles in the system appropriately. The liquid fraction 
must be relatively large so that changes in nucleus size 
significantly affect the density of the surrounding vapor 
and the nucleus cannot evaporate completely because the 
vapor density increases. Eventually, an equilibrium is es- 
tablished, where the nucleus contains n particles while 
the vapor reaches a supersaturated pressure p > p^. 

Farkas [3] pointed out that for a system with N parti- 
cles composed of a nucleus containing n and a supersat- 
urated vapor containing N — n particles, an equilibrium 
between the nucleus and the supersaturated vapor cor- 
responds to the condition ri = n*{p,T). This is due to 
the fact that by definition, the Gibbs energy of nucleus 
formation is maximal for n* {p, T) and thus 
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holds, which implies that for n = n*{p,T), growth and 
decay are equally probable. Since in a nucleation process 
AG„ has a single maximum 0], this equilibrium condi- 
tion uniquely identifies the size of the critical nucleus. 
By minimizing the Hclmholtz energy of the system in an 
NVT simulation, an equilibrium that characterizes the 
maximum of its Gibbs energy is established. The criti- 
cal nucleus model of Reguera-Reiss nucleation theory 27 1 
reproduces these considerations in an explicit form. 

As suggested by Lovett [2^, the fact that a critical 
nucleus 'can only be in (stable) equilibrium with a su- 
persaturated vapour in a system with a finite (small) 
volume' makes these small systems where 'the thermody- 
namic analysis is straightforward and the configurations 
are easily simulated' an attractive topic for molecular 
simulation. Such an approach leads to more accurate 
data on the critical nucleus, e.g. its size n* or surface 
tension 7*, than the usual method of observing growth 
and decay of nuclei in nonequilibrium simulations [l3| . 
because it permits straightforward sampling over a large 
number of time steps. It is also computationally efficient 
since only small systems are considered. The molecu- 
lar simulation of such equilibria is not a novelty in itself 
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Figure 1: (Color) Surface tension of the LJTS fluid over nu- 
cleus size from indirect simulations (A this work, o from pre- 
vious work 18]) and following standard CNT (dashed lines), 
the LFK modification of CNT (dots) , and the new SPC mod- 
ification of CNT (solid lines). The LFK model, which de- 
pends on the supersaturated vapor pressure, weis evaluated 
at 5 = 2.86, 2.28, 1.62, and 1.17 for T = 0.65, 0.70, 0.80, and 
0.95 e/ks, respectively. 



[isl [29l . [sot ' tiut no implications for critical nuclei were 
drawn from these studies in the past. However, Talan- 
quer 



3l| used a similar approach based on density func- 



tional theory for calculating the free energy of formation 
and the interfacial density profile of critical nuclei. 

Simulations in the canonical ensemble based on these 
considerations can contribute to the study of nucleation 
processes indirectly, by reproducing vapor-liquid equilib- 
ria instead of the condensation itself. Such indirect sim- 
ulations were conducted for small systems (total number 
of particles N < 2 x 10'') and properties of the criti- 
cal nucleus at moderate supersaturations were obtained, 
complementing data from previous work . Nuclei with 
10^ < n < 10^ particles were inserted into saturated or 
moderately supersaturated vapor phases. The nucleus 
size was tracked by applying a version of the cluster cri- 
terion of Rein ten Wolde and Frenkel [s^l where a particle 
is considered as belonging to the nucleus if it has at least 
four neighbors within a radius of r < 1.5cr. 

Surface tension and size of the critical nucleus were 
determined for six temperature values between 0.65 and 
0.95 e/fc_B, cf. Fig. [T] as well as Tabs. HI and HIl As usual, 
all numerical results are given in terms of the mass m of a 
single particle and the two potential parameters a (size) 
and s (energy). The planar interface surface tension 700 
of the LJTS fluid, given by [S] 



2.08 1 



T 



(13) 



with Tc = 1.0779 e/kB, is represented by horizontal lines 
in Fig. [1] because standard CNT assumes the surface ten- 
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Figure 2: (Color) Critical nucleus size of the LJTS fluid 
over supersaturated pressure from indirect simulations (A 
this work, o from previous work and following standard 
CNT (dashed lines), the LFK modification of CNT (dots), 
and the new SPC modification of CNT (solid lines). 



sion of the curved interface to be size independent and 
therefore equal to 700 ■ The surface tension, calculated 
from the normal component of the Irving-Kirkwood pres- 
sure tensor [s^l, is significantly reduced for small nuclei 
when compared to 700- Standard CNT neglects this ef- 
fect, and as shown in Fig. [1] the LFK expression for C„ 
given by Eq. (|lip may even lead to values of 7„ for small 
nuclei which are unphysical, i.e. 7„ > 700, and increase 
for smaller nuclei. 

Simulation results from our group for n* are compared 
in Fig. [2] to theoretical values. The LFK modification 
is in better agreement with simulation data over the en- 
tire studied temperature range than standard CNT which 
consistently underestimates n* and leads to particularly 
large deviations at high temperatures. 



DIRECT APPROACH: MD SIMULATION OF 
HOMOGENEOUS NUCLEATION 

A series of direct simulations of the nucleation pro- 
cess was conducted in the canonical ensemble using the 
program Isi [11]. The system size was relatively large 
(10^ < N < 10^) and a hybrid cluster criterion was used 
to detect the nuclei. This criterion combines geomet- 
ric and energetic approaches with a connectivity anal- 
ysis based on graph theory [3^. The nucleation rate 
was determined by defining a threshold i and counting 



the number Ji of nuclei containing at least i particles 
that emerge per volume and time [lOj . According to this 
method, seven nucleation rate isotherms were obtained 
for temperatures between 0.65 and 1 e/ks- 

Nucleation rates Ji for different threshold values i are 
compared to theoretical predictions in Fig. [3] as well as 
Tabs, mil and IIVI The values of Ji are only valid ap- 
proximations of the actual nucleation rate if they remain 
roughly constant for increasing i [lO]. This is the case for 
all temperatures except 0.95 and 1 e/ks where n* is prob- 
ably larger than all of the chosen threshold values. With 
n* ^ i, the rate of formation for nuclei with i or more 
particles does not correspond to the nucleation rate, but 
rather to the velocity at which the metastable equilibrium 
which precedes nucleation is established. For instance, 
at 0.95 e/fcs and a supersaturation of 1.226, standard 
CNT predicts a critical nucleus with 173 particles (LFK: 
n* = 293); thus, the value of Jioo = 3 x em,-^ / a'^ 



Table I; Size of the critical nucleus (in number of particles) 
and its surface tension (in units of e/a^) at low temperatures 
(given in units of e/ks) from simulation in comparison to 
theories - bold values are taken from earlier work qg] 



T 


5 


n* 


7* 


n*(CNT) 


n*(LFK) 


n*{SPC) 


0.65 


1.226 


10500 


0.630 


7800 


7800 


8500 


0.65 


1.337 


3400 


0.610 


2700 


2700 


3000 


0.65 


1.420 


1700 


0.585 


1600 


1500 


1700 


0.65 


1.461 


1300 


0.590 


1200 


1200 


1400 


0.65 


1.512 


1000 


0.578 


950 


940 


1100 


0.65 


1.594 


690 


0.553 


670 


660 


750 


0.65 


1.599 


720 


0.563 


650 


640 


730 


0.65 


1.813 


340 


0.501 


320 


320 


360 


0.65 


1.856 


250 


0.527 


290 


280 


320 


0.70 


1.179 


15600 


0.535 


8900 


8700 


9600 


0.70 


1.306 


2300 


0.507 


2100 


2100 


2300 


0.70 


1.420 


1000 


0.474 


930 


930 


1100 


0.70 


1.474 


740 


0.463 


690 


690 


790 


0.70 


1.586 


420 


0.418 


410 


420 


470 


0.70 


1.621 


430 


0.408 


360 


370 


410 


0.70 


1.722 


230 


0.363 


260 


260 


290 


0.70 


1.816 


130 


0.319 


200 


190 


220 


0.70 


1.869 


150 


0.309 


170 


170 


190 


0.80 


1.125 


9700 


0.355 


8000 


8300 


9100 


0.80 


1.179 


3900 


0.334 


2900 


3100 


3400 


0.80 


1.227 


1600 


0.325 


1500 


1600 


1800 


0.80 


1.264 


1100 


0.308 


1000 


1100 


1200 


0.80 


1.301 


750 


0.264 


730 


790 


870 


0.80 


1.352 


610 


0.270 


490 


530 


580 


0.80 


1.366 


500 


0.251 


450 


490 


520 


0.80 


1.460 


200 


0.190 


260 


280 


290 


0.80 


1.518 


250 


0.152 


190 


220 


220 
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Table II: Size of the critical nucleus (in number of particles) 
and its surface tension (in units of e/cr'^) at high temperatures 
(given in units of e/kB) from simulation in comparison to 
theories - bold values are taken from earlier work [1^ 



T 


S 


n* 


7* 


n 


"(CNT) 


n*(LFK) 


71*(SPC) 


0.85 


1.102 


8800 


0.268 




7100 


7700 


8400 


0.85 


1.136 


3000 


0.250 




3200 


3500 


3800 


0.85 


1.152 


2200 


0.222 




2400 


2600 


2800 


0.85 


1.168 


2200 


0.240 




1800 


2000 


2200 


0.85 


1.227 


840 


0.202 




800 


920 


970 


0.85 


1.253 


1300 


0.224 




600 


700 


730 


0.85 


1.264 


680 


0.157 




540 


630 


650 


0.85 


1.340 


450 


0.127 




290 


340 


340 


0.85 


1.424 


250 


0.072 




170 


210 


190 


0.90 


1.086 


7000 


0.191 




5600 


6400 


6800 


0.90 


1.111 


3700 


0.173 




2700 


3200 


3400 


0.90 


1.134 


2100 


0.165 




1600 


2000 


2000 


0.90 


1.182 


1400 


0.132 




700 


890 


880 


0.90 


1.198 


930 


0.118 




570 


730 


700 


0.90 


1.201 


1000 


0.139 




550 


700 


680 


0.90 


1.223 


650 


0.033 




410 


540 


511 


0.95 


1.067 


7200 


0.114 




4400 


5600 


5700 


0.95 


1.074 


6100 


0.107 




3300 


4300 


4300 


0.95 


1.077 


5400 


0.102 




3000 


3900 


3900 


0.95 


1.082 


3400 


0.092 




2400 


3300 


3200 


0.95 


1.086 


4800 


0.100 




2200 


2900 


2900 


0.95 


1.099 


2600 


0.084 




1400 


2000 


1900 


0.95 


1.104 


1900 


0.059 




1300 


1800 


1700 



obtained under these conditions, cf. Tab. llVi does not 
describe nucleation but equilibration. Some other results 
describe the transition between both regimes, such as 
J400 = 1 X IQ-'^VenF^/a^ at T = 1 e/fcs and S = 1.106, 
where the critical size according to CNT is 328 (LFK: 
n* = 755). 

In those cases where Ji clearly represents the actual 
nucleation process, consistent deviations were found for 
standard CNT which underestimates J by two orders of 
magnitude in all cases. Such an accuracy should be inter- 
preted as a confirmation of standard CNT for the LJTS 
fluid. The LFK modification is in better agreement with 
simulation data at low temperatures but leads to larger 
deviations of J at high temperatures. 



SURFACE PROPERTY CORRECTED 
MODIFICATION OF CNT 



As the preceding sections show, standard CNT only 
predicts the nucleation rate of the LJTS fluid with an ac- 
ceptable accuracy, leading to deviations for n*\ the LFK 
modification provides excellent predictions for the critical 



nucleus size but not for the temperature dependence of J. 
Both theories assume an inappropriate curvature depen- 
dence of the surface tension, although for this essential 
property of inhomogeneous systems a qualitatively cor- 
rect expression is known since the 1940s [2^. With the 
collected simulation data on 7*, n*, and J over a broad 
range of temperatures, enough quantitative information 
is available to formulate a more adequate modification of 
CNT. 

To correlate the simulation results for 7*, Eq. PH)) as 
proposed by Tolman [23[ was chosen. The quotient 5/R 
was assumed to scale with n"^/'^, and a fit to the data 
shown in Fig. [T]as well as Tabs. HI and UIl vields 



0.7 



1 - T/Tc 



-0.9 



1/3 



(14) 



It can be seen in Fig. [T] that the two-parameter fit given 
by Eqs. pUj) and is sufficient to reproduce both tem- 
perature and size dependence of the surface tension. 

The new SPC modification of CNT is based on a size 
dependent term for the surface tension, given by Eqs. (jlOp 
4]) for the LJTS fluid. From theoretical considera- 
and from simulation 36'] it is further known that 
the number of particles n in the nucleus is insufficient as 
a reaction coordinate for nucleation. In particular, the 
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Figure 3: (Color) Nucleation rate of the LJTS fluid over 
supersaturated pressure from the present direct simulations 
for different threshold values (• i = 25, o i — 50, □ i £ 
{75, 100}, A i > 150) and following standard CNT (dashed 
lines), the LFK modification of CNT (dots), and the new SPC 
modification of CNT (solid lines). 
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Table III: Nucleation rate (in units 
temperatures (given in units of elks) 
lations in comparison to theories 



of ^ era ^/o"^) at low 
from the present simu- 
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external shape and ~ in case of liquid-solid nucleation - 
the internal structure of the nuclei must be taken into 
account [s^ [s^l- An MD based analysis of the emerg- 
ing crystals in a supercooled LJ liquid by Trudu et al. 
14| showed that for liquid-solid nucleation, the nuclei 
are significantly anisotropic: the longest axis of the ellip- 
soids used to approximate the crystal surface was found 
to be about 1.5 times longer than the shortest axis. 

The nonsphericity of the nuclei in a supersaturated 
vapor, due to fluctuations of the phase boundary, is rep- 
resented by a size dependent steric coefficient s„ with 
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in the present SPC modification of CNT. The tempera- 
ture and size dependence of s„ was accounted for by a 
two-parameter fit 



Table IV: Nucleation rate (in units of Vem^ /a"^) at high 
temperatures (given in units of e/kB) from the present simu- 
lations in comparison to theories 
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adjusted to all simulation results for n* and J. This cor- 
responds to an effective increase of the radius according 
to capillarity theory i?cap 



0.85(l-T/Tc)-i + (n/75) 
1 + (n/75)i/3 
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^3/(47rp^) by a factor of 
As Fig. [3] shows, this increase is similar in magni- 
tude to phyiscal properties that express the size of the 
phase boundary, such as the Tolman length 5 and the in- 
terface thickness on the vapor side determined from 
the average density profile of the nucleus. Both prop- 
erties were studied for the LJTS fluid in previous work 

For all temperatures above 0.162 e/kB, which is far 
below the triple point temperature (about 0.61 e/fcs), 
Eq. HH) leads to s„ > 1. For T ^ T^, both the steric 
coefficient and the thickness of the vapor-liquid interface 
diverge. It should be noted that in the case of fluids 
that cannot be accurately modeled by the truncated and 
shifted LJ potential, the surface tension of small nuclei 
and the thickness of the phase boundary may require a 
different set of parameters for Eqs. and 

As can be seen from Figs.[T]~|31 the SPC modification 
consistently covers the LJTS simulation results for 7*, n* , 
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and J. At the transition between nucleation and spin- 
odal decomposition, the nucleation rates from simulation 
exceed the prediction by about an order of magnitude, 
which may be due to the particular energy landscape of 
such processes [sl] . The apparent difference between the 
values of Ji and all theories at temperatures of 0.95 and 
1 e/ks is due to the fact that as discussed above, almost 
all of these values describe the velocity of equilibration 
instead of nucleation. 



APPLICATION TO ARGON 

Fluid argon can be represented accurately by an LJTS 
molecular model, the corresponding potential parameters 
were determined as cr = 3.3916 A and e/ks — 137.90 K 
in previous work Conveniently, experimental nu- 

cleation data are available, however, usually at very low 
temperatures below the triple point. The onset pressure 
Pom defined as the pressure where the nucleation rate ex- 
ceeds a certain minimal value Jon, was determined for the 
homogeneous nucleation of argon by Pierce et al. [s^ , Za- 
horanski et al. 0, EH, and Hand et al. The onset 

nucleation rate Jon depends on the experimental setup 
and was provided (or estimated) by the authors of the 
respective studies, except for Pierce et al. f39l . where we 



assumed Jon — lO'^^ m^'^s^^ as given by Hand [8| for the 
nucleation rate detected by supersonic nozzles. 

The comparison between theory and experiment is in- 
conclusive and appears contradictory, cf. Tab.lVl Results 
obtained by Pierce et al. [s^ tend to confirm CNT. Za- 
horanski et al. [40| observed a nucleation onset at much 
lower pressures than all theories. A second study by Za- 
horanski et al. [Jli] agrees best with the SPC modifica- 
tion, whereas the correlation proposed by Hand et al. 
confirms LFK. However, it should be pointed out that 




T = 0.90 s/L 



4000 8000 12000 

nucleus size in number of particles 

Figure 4: (Color) Steric radius factor ^/s^ in dependence 
of the nucleus size n at temperatures of 0.75 and 0.90 e/fcs 
(solid lines); values from previous work [l^ are shown for 1 
-|- 5 /Reap (triangle up) and 1 -I- D'^/Rca.p (triangle down) at 
0.75 (filled symbols) and 0.90 (empty symbols) e/fcs 



Table V: Nucleation onset pressure pon (in units of kPa) for 
argon at low temperatures (in units of K) from experimental 
data in comparison to the pressure where the assumed onset 
nucleation rate Jon (in units of m~^s~^) is reached according 



to theories; the data of Pierce et al. [3E 
graphical form only 



were published in 
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at such very low temperatures (60 to 70% of the triple 
point temperature) desublimation processes take place. 
Therefore, such severe extrapolations to other parts of 
the phase diagram are questionable. 



CONCLUSION 



The present simulation results show for the LJTS fluid 
that standard CNT underpredicts the nucleation rate J 
at high supersaturations by about two orders of magni- 
tude. The critical nucleus size n* is also lower by up to 
a factor of two at moderate supersaturations. The LFK 
modification is an improvement with respect to n*, but 
leads to large deviations for J at high temperatures. A 
surface property corrected modification of CNT was pre- 
sented that takes into account the lower surface tension 
of small nuclei and their nonsphericity, effects ignored 
by standard CNT. This modification consistently repro- 
duces simulation data for 7*, n* , and J over a wide range 
of states. 
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